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1 ABSTRACT 

o 

We report the discovery of type-I X-ray bursts from the low-mass X-ray binary 
4U 1708-40 during the 100 ks observation performed by BeppoSAX on 1999 August 
<i ■ 15-16. Six X-ray bursts have been observed. The unabsorbed 2-10 keV fluxes of the 

Oh bursts range from ~ (3 — 9) xlO -10 erg cm _2 s -1 . A correlation between peak flux 

and fluence of the bursts is found, in agreement with the behaviour observed in other 

■ similar sources. There is a trend of the burst flux to increase with the time interval 
^ | from the previous burst. From the value of the persistent flux we infer a mass accre- 
03 ■ tion rate Af~7x 10 -11 M Q /yr, that may correspond to the mixed hydrogen/helium 

burning regime triggered by thermally unstable hydrogen. We have also analysed a 
BeppoSAX observation performed on 2001 August 22 and previous RXTE observa- 

■ tions of 4U 1708-40, where no bursts have been observed; we found persistent fluxes of 
5h ' more than a factor of 7 higher than the persistent flux observed during the BeppoSAX 

observation showing X-ray bursts. 

Key words: accretion, accretion discs - stars: individual(4U 1708-408) - stars: 
neutron - X-rays: general 



1 INTRODUCTION 

Many low-mass X-ray binaries (hereafter LMXBs) are 
known to show type-I X-ray bursts which have proven to 
be important diagnostic tools for these systems (see Lewin, 
van Paradijs & Taam 1995 for a review). X-ray bursts are 
thought to originate from thermonuclear flashes, caused by 
thermal instabilities, in the accreted matter on the surface 
of weakly magnetized neutron stars (NS). The time interval 
between the bursts ranges from tens of minutes to days. In 
a source the burst activity can stop for periods from days to 
months (but also years, as, for example, in the case of tran- 
sient sources). In some cases, a relation has been observed 
between burst properties and persistent flux. As the persis- 
tent flux increases the recurrence time between bursts tends 
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to increase (typically from a few hours to more than a day) 
and the bursts become less energetic (e.g., Hoffman, Lewin 
& Doty 1977). For instance, in the case of 4U 1705-44 the 
persistent flux increased by a factor ~ 2 while the burst in- 
tervals increased by a factor of ~ 4 (Langmeier et al. 1987; 
Gottwald et al. 1989). When the persistent flux increases 
further the bursts disappear completely. This behaviour is 
observed in several sources like, for example, MXB 1659- 
29 (Lewin, Hoffman & Doty 1976), GX 3+1 (Makishima et 
al. 1983) and EXO 0748-673 (Gottwald et al. 1986). How- 
ever, in other sources (e.g. 4U 1820-30: Clark et al. 1977; 
4U 1728-34: Basinska et al. 1984) the opposite behaviour 
was sometimes observed: the recurrence time between bursts 
decreased while the persistent flux increased. Nevertheless, 
at least in the case of 4U 1820-30 (Clark et al. 1977), the 
bursts (again) disappear for higher acretion rates, when the 
source is in its high state. 
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Current theories of type-I X-ray bursts (Fujimoto et al. 
1981; Fushiki and Lamb 1987; see Bildsten 2000) predict 
three different regimes in mass accretion rate for unstable 
burning: i) mixed H/He burning triggered by thermally un- 
stable H ignition, at low accretion rates (M < 2 x 10 -10 
M a /yr) ; ii) pure He shell ignition following steady H burn- 
ing, at intermediate accretion rates (2 x 10~ 10 M Q /yr < 
M < 4 - 11 x 10~ 10 M /yr); iii) mixed H/He burning trig- 
gered by thermally unstable He ignition at high accretion 
rates (4 - 11 x 10~ 10 M /yr < M < 2 x 10" 8 M /yr). At 
M even higher steady He burning occurs, the nuclear fuel 
for the burst is depleted, causing the bursts to disappear 
completely. On the basis of these theories, an anticorrela- 
tion between persistent flux and burst recurrence time is 
expected; since as M increases a lower amount of time is 
needed to accumulate the critical amount of fuel necessary 
to begin the burst, the recurrence time between bursts is 
expected to decrease with increasing M. This is contrary 
to what is observed (e.g., Hoffman, Lewin & Doty 1977; 
van Paradijs, Penninx & Lewin 1988). There are basically 
two interpretations for these observations, van Paradijs et 
al. (1988) found evidence of a continuous stable burning of 
accreted nuclear fuel in burst sources; the observed correla- 
tion between persistent emission and burst recurrence time 
can be explained as the rate of this stable burning, increas- 
ing, removes nuclear fuel for thermonuclear flash. Recently, 
Bildsten (2000) proposed that another solution of this mis- 
match between theory and observations can be obtained by 
relaxing the hypothesis of spherical accretion and assuming 
that the accretion of fresh material only occurs on a fraction 
(e.g., the equator) of the NS surface. The crucial parameters, 
in this case, is the accretion rate per unit area: if the covered 
area increases enough with increasing M, the accretion rate 
per unit area indeed decreases. Indeed, the whole scenario is 
not completely clear yet, since recent results (e.g. Muno et 
al. 2000; Franco 2001; van Straaten et al. 2001) suggest that 
there is not a unique trend between M and burst properties 
for all the sources. 



Type-I X-ray burst profiles depend strongly on photon 
energy: decays are much shorter at high photon energies 
than at low photon energies (Lewin et al. 1995). This soften- 
ing of the burst spectrum during the decay results from the 
cooling of the NS photosphere. If the luminosity during the 
burst reaches the Eddington limit LEdd, the NS photosphere 
expands. Since for a blackbody (see below) Lb oc R 2 T^f, 
when the radius of the photosphere R expands, the effec- 
tive temperature T e / / decreases (Tawara et al. 1984; Lewin, 
Vacca & Basinska 1984; Vacca, Lewin & van Paradijs 1986). 



In this paper we report on the discovery of type-I X-ray 
bursts from the LMXB 4U 1708-40 (Forman et al. 1978) 
during a BeppoSAX observation. Not much is known about 
this system. A first analysis of the 2-12 keV X-ray spectrum 
of 4U 1708-40 was made by Warwick et al. (1988) using 
EXOSAT observations during which the 2-6 keV source flux 
was ~ 10 -9 erg s _1 cm~ 2 . The spectrum was well fitted by a 
power law with a spectral index T = 2.2 ± 0.2 and a column 
density N H ~ 3 x 10 22 cm -2 . 



Table 1. Best-fit parameters and la errors of the spectra of the 
persistent emission during the 1999 August 15-16 (SAX 1999), 
2001 August 22 (SAX 2001) BeppoSAX observations (0.12- 
10 keV) and the 2000 June 18 PCA/RXTE observation (2.5- 
25 keV). The 1999 BeppoSAX data are fitted with an absorbed 
power law or a CompTT model plus a Gaussian emission line, the 
2001 BeppoSAX data with an absorbed power law plus a black- 
body component and the RXTE data with an absorbed cutoff 
power law plus a blackbody. Nh is the equivalent hydrogen col- 
umn density in units of 10 22 cm -2 , N p i is the normalisation of 
the power law in units of 10~ 2 photons/keV/cm 2 /s at 1 keV, T is 
the photon index of the power law, To is the temperature of the 
soft seed photons for the Comptonisation, T e is the electron tem- 
perature, t is the optical depth of the spherical scattering cloud, 
N T t is the CompTT normalisation (xl0~ 2 ) in XSPEC vll.1.0 
units (see Titarchuk 1994a, b), N cu t is the normalisation of the 
cutoff power law in units of 10 — 1 photons/keV/cm 2 /s at 1 keV, 
E cu t is the cutoff energy, N bb is the normalisation of the black- 
body in units of L37/TJ 2 (where L37 is the source luminosity in 
units of 10 37 erg/s and D10 is the distance of the source in units 
of 10 kpc) and T bb the color temperature of the blackbody. Eij ne 
is the centroid energy, cr; ine is the width and eqw is the equiva- 
lent width of the Gaussian emission line. Ftot is the unabsorbed 
2-10 keV flux in erg s _1 cm -2 . 



SAX 


Parameters 


Power law 


CompTT 


1999 




+ Gaussian 


+ Gaussian 




N H 


2.93 ±0.08 


2.45 ±0.14 




r 


2.42 ±0.02 






N p i (LECS) 


8.54 ±0.34 






N p i (MECS) 


10.72 ±0.31 






kT (keV) 




0.32 ±0.03 




kT e (keV) 




5.30 ±0.9 




T 




3.3 ±0.3 




N TT (LECS) 




1.4 ±0.2 




N TT (MECS) 




1.7 ±0.2 




E line (keV) 


6.5 ± 0.1 


6.5 ±0.1 




a Hne (keV) 


0.9 


0.9 




eqw (eV) 


344 


286 




Ftot 


1.2 x lO" 10 


1.2 x lO" 10 


SAX 


Parameters 


Power law 




2001 




+ blackbody 






N H 


2.9 (fixed) 






r 


2.68 ±0.29 






N p i (MECS) 


41.88 ±9.09 






kT bb (keV) 


1.31 ±0.03 






N bb (MECS) 


1.02 ±0.11 






F tot 


8.9 x lO" 10 




RXTE 


Parameters 


Cutoffpl 








+ blackbody 





N H 

r 

N cut 

Ecut (keV) 
kT bb (keV) 
N bb 
Ftot 



2.9 (fixed) 
1 44+013 

3.30 ± 0.38 
4.84 ±0.36 
1.27 ±0.02 
0.36 ±0.03 
7.4 x 10- 10 



Discovery of type-I X-ray bursts from the LMXB 4U 1708-4 3 




10 20 30 40 50 

Time (sec) 



^ - 




10 20 30 40 50 

Time (see) 



Figure 1. Upper panel: light curve with 1 s time resolution of the 
first burst. Lower panel: hardness ratio (HR= 4 — fO keV/2 — 4 
keV) curve of the first burst with 2 s resolution. The shape of the 
burst, the hardening during the rise and the softening during the 
decay are typical of a type-I X-ray burst. 



and MECS combined to analyse the persistent emission, the 
second, third, fourth and fifth bursts, and only the MECS 
to analyse the first and the sixth bursts, since the LECS 
was off in these time intervals. The BeppoSAX observation 
is periodically interrupted due to earth occultations and the 
passage through the South Atlantic Anomaly. These 'gaps' 
last on average ~ 2 ks after ~ 4 ks of on-source observation. 
The BeppoSAX instruments were off more often in the last 
part of the observation (namely between the fourth and the 
sixth burst) so that the on-source exposure time is ~ 50 ks. 

We have also analysed a 150 ks BeppoSAX observation 
of 4U 1708-40 performed on 2001 August 22 (with an on- 
source exposure time of ~ 73 ks) and 14 (non contiguous) 
observations from the public RXTE data archive, five per- 
formed in 1997 and nine in 2000, with a total on-source ex- 
posure time of ~ 60 ks. For the RXTE observations we have 
used data taken with the Proportional Counter Array (PCA, 
which consists of five co-aligned Proportional Counter Units, 
PCUs, sensitive in the energy range 2-60 keV; Zhang et al. 
1993). Starting from 2000 May 12, the propane layer on 
Proportional Counter Unit (PCU0), which functions as an 
anticoincidence shield for charged particles, was lost. There- 
fore we excluded the PCU0 data from the 2000 June spec- 
trum we have analysed. The PCA observations background, 
estimated using peabekest v2.1e was also subtracted. 

We produced light curves of the 2001 BeppoSAX obser- 
vation and of all the RXTE observations; in all these data 
we did not find any burst. Although a small variation in av- 
erage count rate (around 20%) occurred in the RXTE data 
between the 1997 observations and the 2000 observations, 
the light curve of each observation shows steady persistent 
emission without any significant variation in count rate; we 
analysed the spectrum of just one of the RXTE observa- 
tions available, as a representative case. To analyse the PCA 
spectrum we have used Standard2 data in the energy range 
2.5-25 keV and produced the detector response matrix with 
pcarsp v7.11. A systematic error of 1% was added to ac- 
count for residual uncertainties in the detector calibration. 
We have used XSPEC vl 1.1.0 to fit the spectra. 



2 OBSERVATIONS AND DATA ANALYSIS 

We have analysed a BeppoSAX observation of 4U 1708- 
40, performed on 1999 August 15-16 for a total observa- 
tion time of ~ 100 ks. We report on results from the Low 
Energy Concentrator Spectrometer (LECS, a thin- window 
position-sensitive gas scintillator proportional counter with 
extended low energy response, 0.1-10 keV, and a field of 
view, FOV, of 20'; Parmar et al. 1997) and the Medium 
Energy Concentrator Spectrometers (MECS, position sensi- 
tive gas scintillator proportional counters operating in the 
1.3-10.5 keV band, with a FOV of 30'; Boella et al. 1997b). 
The source was not detected in the Phoswich Detection Sys- 
tem (PDS, energy range 13-200 keV; Frontera et al. 1997). 
For both imaging instruments we selected the data in cir- 
cular regions of radius 8' centered on the source. We used 
the standard response matrices and background files (1997 
version for MECS and 2000 for LECS) for the spectral anal- 
ysis. To analyse the BeppoSAX spectra of the source we 
selected the energy range 0.12-4 keV for LECS and 1.8-10 
keV for MECS. In the 1999 August 15-16 BeppoSAX obser- 
vation six X-ray bursts are detected. We have used LECS 



3 RESULTS 

In the 1 s time resolution LECS and MECS light curves of 
the 1999 observation we see persistent emission, at a count 
rate of ~ 2 - 4 counts/s in the MECS light curve. The light 
curve also shows six X-ray bursts (the times at which each 
burst occurred are reported in Table 2). The bursts last 
~ 40 s. Excluding the bursts, the persistent emission does 
not show significant intensity variations during the whole ob- 
servation. Five bursts out of six exhibit a rapid rise ( < 5 sec) 
followed by a slower decay, whereas one of the bursts (the 
fourth) has a triangular shape with approximatly equal rise 
and decay times of 20 sec. In Fig. 1 we show the light curve 
in the MECS range (upper panel) and the corresponding 
hardness ratio (HR= 4—10 keV/2 — 4 keV, bottom panel) 
of the first burst, the one observed with the best statis- 
tics (the behaviour of the other bursts are comparable with 
this). The hardening during the rise and the softening dur- 
ing the decay, together with the shape of the bursts, identify 
4U 1708-40 as a type-I X-ray burster. 
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Figure 2. Upper panel: spectrum from BcppoSAX data (LECS 
and MECS) of the persistent emission of 4U 1708-408 in the 
range 0.12—10 keV, averaged over more than 45 ks of observation, 
and the best-fit model (absorbed power law) reported in Table 1. 
Lower panel: residuals in units of a with respect this model. 
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Figure 3. Upper panel: RXTE/PCA spectrum of the persistent 
emission of 4U 1708-408 in the 2.5-25 keV energy band, with 
the best-fit model (absorbed cutoff power law and blackbody, 
Table 1; residuals around 6-7 keV are still present: see § 3.1). 
The spectrum is averaged over the whole ~ 3 ks observation of 
2000 June 18. Lower panel: residuals in units of cr with respect 
this model. 

3.1 Persistent Emission 

We produced LECS and MECS spectra of the persistent 
emission of the 1999 August 15-16 observation of 4U 1708- 
40 excluding intervals of 200 s around (starting about 70 s 
before) each of the six detected bursts. The spectrum in 
the range 0.12-10 keV is shown in Fig. 2. We fit the spec- 
trum with a power law corrected for photoelectric absorp- 
tion and a Gaussian emission line at 6.5 ± 0.1 keV which 
gives a reduced x 2 of 1.2 for 420 d.o.f. (the reduced \ 2 
without the Gaussian emission line is 1.3 for 423 d.o.f.). 
We find a high equivalent hydrogen column density, Nh = 
(2.93T0.08) x 10 22 cm~ 2 , consistent with the position of the 
source in the direction of the galactic center, and a power 
law photon index V = 2.42 ± 0.02. No significant thermal 
component (i.e., blackbody) is found (with an upper limit 
on the flux of ~ 2 x 10~ 12 erg cm _2 s _1 , fixing the black- 
body temperature to 1.3 keV, see below). The 2-10 keV 
flux is 1.2 x 10~ 10 erg cm~ 2 s -1 , corresponding to a lumi- 



nosity of ~ 9 x 10 35 erg s _1 at a distance of 8 kpc (see Sec- 
tion 4). We obtain a fit of similar quality (reduced \ 2 — 1 -2 
with 421 d.o.f.) using the Comptonisation model CompTT 
(Titarchuk 1994a, b), with a slightly lower column density, 
Nh = (2.56 ± 0.14) x 10 22 cm -2 , and a Gaussian emission 
line at 6. 5 ±0.1 keV. In Table 1 we show the fit parameters of 
the persistent emission corresponding to these two models. 

For comparison we have also analysed the 4U 1708-40 
spectrum during one of the RXTE observations (the 3 ks 
observation performed in 2000 June 18, obsID 50060-04-01- 
00). A simple absorbed (the column density was fixed to 
Nh = 2.9 x 10 22 cm~ 2 , in accordance with the BeppoSAX 
spectrum) power law is not able to fit the spectrum in the 
whole energy range. A blackbody component with a color 
temperature of ~ 1.3 keV is needed to adequatly fit the data 
(an F-test gives a probability of chance improvement of the 
fit for the addition of this component of 2 x 10~ 7 ). A high 
energy cutoff at relatively low energy (E cut ~ 5 keV) is also 
necessary (reduced \ 2 = 0.9 with 44 d.o.f.). We show this 
spectrum, together with residuals in units of o with respect 
to the best-fit model, in Fig. 3 and the best-fit parameters 
in Table 1. We measure a 2-10 keV unabsorbed flux of ~ 
7.4 x 10~ 10 erg cm^V, a factor ~ 7 times higher than 
the flux of the persistent emission we measure during the 
BeppoSAX observation. In Fig. 3 we note residuals around 
6-7 keV. For comparison with the spectrum during the 1999 
August 15-16 BeppoSAX observation, we tried to add a 
Gaussian emission line in the range 6.4-6.9 keV. We find 
a Gaussian line at ~ 6.75 keV with an equivalent width of 
87 eV. The parameters of the other components does not 
change significantly. 

We have also analysed the MECS spectrum of the 
BeppoSAX observation performed on 2001 August 22. The 
best-fit is obtained using a power law with a slope of 
T = 2.68 ± 0.28, and a blackbody with a color tempera- 
ture of kT= 1.31 ± 0.03 keV (reduced \ 2 = 1-07 with 172 
d.o.f.; Table 1). The spectrum is compatible with a Gaus- 
sian emission line between 6.4 and 6.7 keV, similar to the 
one used to fit the source spectrum during the previous Bep- 
poSAX and RXTE observations, although this component 
is not required in this case. The unabsorbed 2-10 keV flux 
is 8.9 x 10~ 10 erg cm _2 s _1 . This model is consistent with 
the best-fit model of the PCA 2-10 keV spectrum, that is 
well fit with a power law with aT~ 2.77 and a blackbody 
with a color temperature of ~ 1.3 keV. 



3.2 X-ray Bursts 

We have analysed the six X-ray bursts of the 1999 Bep- 
poSAX observation. We do not have enough statistics to 
select different intervals during the bursts and analyse the 
rise and decay spectra separately. Therefore we have anal- 
ysed six spectra (one for each burst) each averaged over the 
whole ~ 40 s burst duration. From these spectra we have 
subtracted the spectrum of the persistent emission and fit- 
ted them with a blackbody component. In Table 2 we show 
the results of the fits for each of the six bursts. Note that, 
because the spectra of the bursts soften significantly dur- 
ing the observation (see Fig. 1) the fit parameters, such as 
blackbody color temperatures, have to be considered aver- 
age estimates over the bursts. 



Discovery of type-I X-ray bursts from the LMXB 4U 1708-4 5 



Table 2. Blackbody color temperature (kTbb in keV), normalisa- 
tion (N;,;,, in units of L37/Dj where L37 is the source luminosity 
in units of 10 37 erg/s and D10 is the distance to the source in 
units of 10 kpc) and 2-10 keV unabsorbed flux F (xl0~ 10 erg 
cm - 2 s — 1 ) with la errors, from the fit of the averaged spectra over 
40 s of the six bursts with a blackbody component. In the first 
column in brackets is also indicated the time (with 1 s time reso- 
lution) from the beginning of the MECS observation (15 August 
1999 at 05:01:10.5) at which the peak of each burst occurs. 



burst # (time) kTbb 



N bb 



1 (21655 s) 2.01 ±0.16 1.52 ±0.17 9.2 ± 0.6 

2 (29160 s) 1.60 ±0.13 0.75 ± 0.08 5.6 ± 0.4 

3 (39620 s) 1.64 ±0.10 1.12 ±0.09 7.7 ±0.5 

4 (52350 s) 1.75 ±0.12 1.49 ±0.14 9.6 ± 0.7 

5 (57498 s) 1.28 ±0.09 0.68 ± 0.06 5.1 ±0.4 

6 (91064 s) 1.40 ±0.15 0.45 ± 0.06 3.8 ± 0.3 



4 DISCUSSION 

We have analysed two BeppoSAX observations (taken in 
1999 and 2001) and RXTE observations (taken in 1997 and 
2000) of 4U 1708-40 and discovered this source to be an 
X-ray burster. This allows us to classify 4U 1708-40 as a 
NS system. In the 1999 BeppoSAX observation we found 
six bursts. Five bursts have a rapid rise ( < 5 s) with a slow 
decay (~ 35 s), and one has a triangular shape, although it 
shows spectral properties similar to those of the other bursts. 
The fluxes of the persistent emission in the RXTE observa- 
tions and in the 2001 BeppoSAX observation are ~ 7 and 
~ 8 times, respectively, higher than the flux of the persis- 
tent emission in the 1999 BeppoSAX observation. While we 
observe six bursts during the 1999 BeppoSAX observation 
we do not observe any burst during the RXTE observations 
and during the 2001 BeppoSAX observation. This would be 
in agreement with the general behaviour of X-ray bursters: 
the bursts disappear above a certain flux value due to steady 
nuclear burning (e.g., van Paradijs et al. 1979; van Paradijs 
et al. 1988; see Lewin et al. 1995 for a review). 

Several sources show a correlation between the peak flux 
and the fiuence of the bursts (e.g. Sztajno et al. 1983; Basin- 
ska et al. 1984; Lewin et al. 1987). Some of them also show 
a saturation in the peak flux at high fluences, which has 
been explained by the fact that the luminosity has reached 
a critical value that can be interpreted as the Eddington 
limit luminosity. One of the best examples is MXB 1728-34 
(Basinska et al. 1984; Di Salvo et al. 2000; Galloway et al. 
2002). This property can be used to infer the distance D to 
the source. In Fig. 4 we plot the flux at the peak F p (cal- 
culated using the count rate at the peak of the burst in the 
0.3 s resolution light curve) versus the fiuence Eb (the aver- 
age bolometric flux times the duration of the burst) for each 
of the six bursts of 4U 1708-40. The plot shows an approx- 
imately linear correlation between F p and Eb- Although we 
do not see a saturation of F v in this correlation (and there- 
fore we cannot derive the Eddington luminosity and the dis- 
tance to the source) we can at least give an upper limit to 
the distance assuming the highest peak flux of the bursts to 
be less than L E dd- Since F p x 47rD 2 < L E dd ~ 2.5 x 10 38 erg 
s _1 (assuming M = 1.4 M Q and correcting for gravitational 
redshift; van Paradijs & McClintock 1994), taking the bolo- 




Figure 4. Flux at the peak (F p ) as a function of the fiuence (Eb) 
of the six bursts. The empty triangle represents the burst with a 
triangular shape. 
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Figure 5. Fiuence Ej, of five bursts (second, third, fourth, fifth: 
open circles; sixth: filled circle) as a function of the time interval 
from the previous burst. 



metric flux at the peak of the first burst (F p = 8.6 x 10~ 9 
erg cm~ 2 s _1 ), we obtain D < 16 kpc, not a stringent con- 
straint. Since the source is in the direction of the galactic 
center we assume 8 kpc as the source distance. At 8 kpc 
saturation would occur at 4 times this flux level. 

Note that, as was pointed out by many authors, some- 
times not all the accreted fuel is burned during the burst 
event, and the presence of a residual fuel can have significant 
implications for the lack of regularity of the burst behaviour 
(Ayasly & Joss 1982; Hanawa & Fujimoto 1984; Woosley & 
Weaver 1985; Lewin et al. 1995). To investigate this point 
we have also calculated the values of a=(GM/R)/.E nuc (an 
observational quantity defined as the ratio of the average 
total energy in the persistent emission to that emitted in 
the burst) which is expected to be ~ 40 for a thermonuclear 
burst in the mixed hydrogen/helium burning regime (e.g., 
Lewin et al. 1995). We measure a ~ 40 for all the bursts, 
but one, starting from the first burst. Given that the accre- 
tion rate is consistent with being constant, this implies a 
correlation between the time interval during which matter 
accretes onto the NS and the energy emitted in the following 
burst. The only exception is the last burst for which, assum- 
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ing that the fifth and the sixth bursts are consecutives, we 
find a value for a that is much larger than 40 (a ~ 250). 
In Fig. 5 we plot the fluence of the bursts as a function of 
the time interval between each burst and the previous one. 
There is a clear trend of the fluence to increase with the time 
interval for all the bursts, except the sixth. This trend is in 
agreement with the behaviour of other sources (Lewin et al. 
1995 for a review) and with the expectation that the longer 
the burst interval is, the larger is the amount of nuclear fuel 
available for the burst. The sixth burst seems not to follow 
this trend. Most probably we have missed some bursts in 
between the fifth and the sixth because of the presence of 
gaps (see §2). Indeed, we find a gap about 3000 s before the 
sixth burst, just where, following the correlation in Fig. 5, 
we would expect to find another burst. 

Based on the observed persistent flux and the distance 
to the source, we can estimate the accretion rate and there- 
fore the burning regime expected in this case. We find 
M ~ 7x 10~ n Mq/yt. This corresponds to the mixed hydro- 
gen/helium burning regime triggered by thermally unstable 
hydrogen (Fujimoto et al. 1981; Fushiki and Lamb 1987; 
see also Bildsten 2000) . This is consistent with the value of 
a ~ 40 we find for almost all the bursts. 
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